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ABSTRACT 

The detection of polarized sources in the WMAP 5-year data is a very difficult task. The maps are 
dominated by instrumental noise and only a handful of sources show up as clear peaks in the Q and U 
maps. Optimal linear filters applied at the position of known bright sources detect with a high level of 
significance a polarized flux P from many more sources, but estimates of P are liable to biases. Using a 
new technique, named the filtered fusion technique, we have detected in polarization, with a significance 
level greater than 99.99% in at least one WMAP channel, 22 objects, 5 of which, however, do not have 
a plausible low radio frequency counterpart and are therefore doubtful. Estimated polarized fluxes 
P < 400 mJy at 23 GHz were found to be severely affected by the Eddington bias. The corresponding 
polarized flux limit for Planck/LFI at 30 GHz, obtained via realistic simulations, is 300 mJy. We 
have also obtained statistical estimates of, or upper limits to the mean polarization degrees of bright 
WMAP sources at 23, 33, 41, and 61 GHz, finding that they are of a few percent. 
Subject headings: filters: filters; point sources: catalogues, identifications, polarization 



1. INTRODUCTION 

Studies of Cosmic Microwave Background (CMB) 
anisotropy are a top scientific priority since they address 
the deepest questions about origin, structure, and equa- 
tion of state of the Universe. Given the sensitivity of 
current detectors, the main constraint on our ability to 
accurately map CMB anisotropics is set by foreground 
emissions. While these signals have contaminated, but 
not dominated, temperature maps, they are a far big- 
ger problem for CMB polarization. Also, fighting this 
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contamination is more difficult because we know much 
less about polarization than we do about total int ensity 
emiss ion. The Task Force on CMB Research (Bock et al.l 
I2006D indeed regards a better characterization of polar- 
ized foregrounds as "a key milestone" in their proposed 
roadmap. 

Extragalactic radio sources are the main CMB contam- 
inant on angula r scales belo w 0.5° at frequencies of up 
to ~ 100 GHz (jTucci et al.l 12005). Observational stud- 
ies of high-frequency polarization are still scanty and 
mostly dealing with source s selec ted at lower frequen- 
cies as shown in'Ri cci et al. (I2004D and in the references 
in Table 3 of .Tucci et all (|2004l ). The bhnd AustraUa 
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Telescope 20 GHz (AT20G) survey of the Southern sky 
includes polarization measurements at 20, 8.6 and 4.8 
GH z: data have bee n publish ed for the bright source s am- 
ple (jMassardi et al . 2008; B urke-Spolaor et al.l[200i ). 

The WMAP survey has yielded the first all-sky to- 
tal intensity and polarization surveys at fre quencies from 
23 to 94 GHz. The analysis of 5-yr data (jWright et all 
|2009() showed that, in general, the WMAP detected point 
sources are not strongly polarized. Only 5 (Fornax A, 
Pictor A, 3C 273, Virgo A, and 3C 279) were found to 
have polarization degrees greater than 4% in two or more 
bands. In this paper we plan to complement and im- 
prove on their analysis in two basic respects. On one 
side, we apply a non-blind approach to source detec- 
tion in polarization, exploiting the knowledge of posi- 
tions of the brightest sources in total intensity, and a 
new d etection technique, called the filtered fusion tech- 
nique (jArgiieso e t al, 2 009.,) . taking into account the real 
beam profiles. On the other side, we check the reliabil- 
ity of our estimates of polarized flux densities by com- 
parison with the very high signal- t o-nois e AT20G BSS 
(Massardi et al. 2008). lRicci et al.l (|2004D . and high fre- 
quency VLA calibrator measurements, and by means of 
numerical simulations. We also present an estimate of 
the mean polarization of sources in a total flux density 
limited sample. 

The polarization is measured by the Stokes 
parameters Q, U and V, and the polar- 



ized intensity is P — yQ^ + t/^ + T^^, see 
IKamionkowski. Kosowskv fc Stebbinj (|1997f ) for further 
details. If we consider linear polarization, = 0, 
which is justified by the fact that extragalactic radio 
sources have very low levels of circular polarisation 
(jHoman. Attridge &: Wardiil2001f) . we have to combine 
Q and U maps in an appropriate way to avoid biasing 
the estimates of P. Furthermore, the polarized signal 
is just a small fraction of the total intensity signal, 
of the order of a few percent, which makes it hard to 
detect. In the case of WMAP polarization maps, only a 
couple of sources can be seen by eye. This situation may 
exacerbate the problems with the Eddington (1913) bias 
(sources are more easily detected if they happen to lie 
on top of positive fluctuations, so that, o n average, their 
fluxes are overestimated), highlighted bv lMassardi et al.) 
(|2009) for WMAP sources close to the detection limit in 
total flux. This issue will be further investigated. 

The outline of this paper is as follows. In §[5] we de- 
scribe our method and the sample of WMAP sources, 
bright in total flux, to which it was applied to estimate 
their polarized flux. In §[3] we discuss the results, also in 
comparison with high signal-to-noise ground based mea- 
surements at similar frequencies. Finally, in §|4l we sum- 
marize our main conclusions. 

2. METHODOLOGY 

2.1. The Input Catalog 

We have carried out a systematic investigation of the 
polarized flux o f the 516 sources detected at > 5ct by 
iMassardi et all (|2009( ) in the 5-yr WMAP temperature 
maps, and listed in the NEWPS_5yr_5s catalogue^. 484 

^ http:/ /max. ifca.unican.es/caniego/NEWPS 



of these sources have a clear identification in low fre- 
quency catalogues (including 27 Ga lactic objects) ; the 5 
objects detected in polarization by IWright et al.l ()2009f ) 
belong to this group. The remaining 32 candidate sources 
do not have plausible counterparts in all-sky low radio 
frequency surveys and may therefore be just high peaks 
in the highly non-Gaussian distribution of the other com- 
ponents present in the maps. If they are all spurious, the 
reliability of the sample is 94% and its completeness is 
of 91% above 1 Jy. 

Since one of our goals is to define a sample as large 
as possible of potential calibrators for CMB polarization 
experiments, we have added 3 extended sources (Cygnus 
A, Taurus A, Gas A), not included in the catalog be- 
cause they lie in very noisy regions, close to the Galactic 
equator, but known to be very bright and significantly 
polarized. The full sample (Input Gatalog) is thus made 
of 519 sources. 

2.2. Filtered Fusion 

Methods to extract astrophysical foregrounds from 
multi-frequency GMB maps frequently exploit the 
prior knowledge of their frequency dependence. This 
approach however does not work well for radio sources 
because of the broad variety of their s pectral proper- 
ties i n the releva,nt fre quency range (jMassardi et al.l 
120081: ISadler et al] l2008f ). On the other hand, with 
few exceptions, extragalactic radio sources look point- 
like when observed with the beams used by GMB 
experiments, and therefore have, in the maps, the 
shape of the effective angular response function of 
the instrument. In the literature one can find several 
methods exploiting this property and using linear 
filters to detect point sources. The standard matched 
filter approach has b e en u se d for years (iNailone 
19921: iVikhhnin eFal] 119951: iMalik fc Subramanian 



1997 



Tegmark fc de Oliveira-Gostal 



1998 



2001 



Sanz. Herranz fc Martf nez-Gqnzale: 
Herranz et aLrT2002t iStewartI 120061 ). More recently, 
a multi-frequency approach based on matched filters has 
been elaborated (Herr anz et al.l 120091 : iHerranz fc Sand 
120081). An approach based on optimal wavelets 
^(IVielva et a l. 2001, 2003 : iBarnard et al.ll200l iSanz et all 
120061 iGonzalez-Nuevo et all I2006D was successfu l ly ap- 
plied to WMAP m aps (jLopez Ganiego et al.l l2007t 
IMassardi et al.l 120091 ) as well as to realistic simu- 
lations of Planck maps (jLopez-Ganiego et al.l 120061 : 
iLeach et al.l |2008() . Moreover, filters based on the 
Neyman-Pearson approach, using the distribution of 
maxima, have been proposed (Lopez-Ganiego et al. 
2005a, b) and a Bayesian approach has been developed 
jHobson fc McLachlan 2003; Fcroz fc Hobson 200i; 
' Carvalho. Rocha. fc Hobsonll2009l) . 

In this work, following Argiieso et al. (2009), we use 
the same matched filter over Q and U images. The 
matched filter is a circularly-symmetric filter, ^'(a;; i?, &), 
such that the filtered map, w{R, b) satisfies the following 
two conditions: (1) (w(i?o,0)) = s(0) = A, i.e. w{R,0) 
is an unbiased estimator of the fiux density of the source; 
(2) the variance of w{R, b) has a minimum on the scale 
i?o, i-G- it is an efficient estimator. In Fourier space the 
matched filter writes as: 



MF — — 



1 T{q) 



aA{q)' 



2tt / dqq 



A(g) 



(1) 
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Fig. 1. — Polarization maps at 23 GHz centered on the position of Fornax A (no. 74 in the NEWPS_5yr_5s catalogue). The upper left- 
and right-hand panels show, respectively, the maps of Q and Q before and after filtering with the matched filter, while the middle panels 

show the analogous and U maps. The lower panels show, for the same source, the maps of P = yCp + fj^ (left) and P = + C/^ 

(right). 



where A((7) is the power spectrum of the background and 
r((7) is Fourier transform of the source profile (equal to 
the beam profile for point sources). 

Since, in this application, each patch is centered on 
the position of a source detected in total intensity, we 
describe the source as 

s{x)=At(x), (2) 

where A is its unknown polarized flux density and t{x) is 
its profile. We assume circular symmetry, so that t{x) = 
t{x), X = \x\. For point sources the profile is equal to 
the beam response function of the detector. The WMAP 
beams are not Gaussian and we use the real symmetrized 



radial beam profiles for the different WMAP channels to 
construct our filters. 

The matched filter gives directly the maximum amplifi- 
cation of the source and yields the best linear estimation 
of the flux. As extensively discussed in the literature, 
it is a very powerful tool to detect point sources, but it 
has to be used with care because its performance may 
degrade rapidly if it is not properly implemented. In 
particular, the power spectrum A(q) of the image has to 
be obtained directly from the data. 

The WMAP team has used a matched filter that op- 
erates on the sphere and takes into account the non- 
Gaussian profile of the beam. It is described by eq. ([T]), 
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replacing the flat limit quantities r(q) and A(q) with 
th eir harmonic equivalents bi and Ci. As pointed out 
bv lLopez-Caniego et al.l ()2006D . the use of the Ce's com- 
puted from all-sky maps to construct a global matched 
filter that operates on the sphere is a reasonable first 
approach, but we think that it can be improved by op- 
erating locally. 

We have followed the sch eme named filtered fusion by 
its authors (jArgiieso et al.ll200 9'). For each object in the 
input catalog, and for each WMAP frequency between 
23 and 61 GHz, we have projected two patches (one for 
Q and one for U), each of 14.65 x 14.65 square degrees, 
centered on the source position. We have left aside the 
94 GHz channel because of the normalization problems 
discussed by Lopez-Caniego et al. (2007) and Gonzalez- 
Nuevo et al. (2008). Each patch is made of 128 x 128 
6.87 arcmin pixels (HcalPix Nsidc = 512, IGorski et al.l 
([2005,)). The projection has been done using the CPACK 
library^. Then, each pair of patches has been filtered us- 
ing a matched filter exploiting the power spectrum de- 
termined within each patch. After filtering, the and 
are added together and the square root of the re- 
sulting image has been calculated. The noise bias can 
be removed by subtracting from and the corre- 
sponding noise contribution and cr? . This correction 

turns out to be negligible. In this way we have obtained a 
map of the polarized intensity P within each patch. This 
approach differs from the usual one, where a P map is 
constructed adding together the unfiltered Q"^ and U"^ 
maps, and taking the square root of the resulting map 
(see Fig.[T]). 

We have then looked for the brightest maximum inside 
a circle centered on the center of each patch and covering 
the area of the WMAP beam. Next, we estimated the 
noise level of the patch and the significance of the possible 
detection. Finally, we constructed catalogues containing 
all sources whose polarized fiux P was detected above a 
chosen significance level. 

The reason why we need to use the significance of the 
detection instead of the usual signal-to- noise ratio is that 
the noise does not obey a Gaussian distribution but a 
Rayleigh distribution, since we are dealing with maps 
that have been squared. The significance was derived 
from the distribution of the values of P for the pixels 
within the patch. An example of such distribution is in 
Fig.l 

2.3. Flux and error estimation 

The polarized flux densities and their errors were esti- 
mated in a way similar to that applied for total intensity 
(Massardi et al. 2009). Point sources appear in the image 
with a proflle identical to the beam profile. For example, 
if the beams were Gaussian, the source flux could be ob- 
tained multiplying the flux in the brightest pixel by the 
ratio between the beam and the pixel area, 2tt{Rs)^ / L^, 

where Rs = FWHM/(2V2 log 2) and Lp is the pixel side. 

In our case the beams are not Gaussian and we need to 
calculate this relationship integrating over the real sym- 
metrized beam profile for each channel. In doing that we 
have to take into account that we work with HEALPix 
pixehzation ()G6rski et al.ll2005D at Nside=512. Although 

^ http:/ /astro. ic.ac.uk/mortlock/cpack/ 




Fig. 2. — Histogram showing the distribution of the polarized 
flux P = \/(Q^~+~U^ (in mK per pixel) obtained from a filtered 
WMAP 23 GHz patch of Q and U centered on the position of 
Fornax A (14.65 square degree patch and 128x128 pixels). This 
histogram has been produced with the values of ~ 13500 pixels, 
excluding the flagged ones (see Fig. (Sjl. The vertical lines corre- 
spond to the values of P exceeding those of 99.0%, 99.9%, and 
99.99% of the pixels. In other words, measured values of P at 
these levels have, respectively, 99.0%, 99.9%, and 99.99% probabil- 
ity of not being due to noise spikes. Note that the polarized flux 
of a source is obtained multiplying the flux in the brig htest pixel 
by the ratio between the beam and the pixel area (see ^ 12.31 1. 

the image is centered at the position of the source, af- 
ter the projection to the fiat patch the object does not 
always lie in the central pixel, but may end up in an ad- 
jacent one. Thus, to estimate its flux we make reference 
not to the intensity in the central pixel but to that of the 
brightest pixel close to the center of the filtered image 
within an area equal to that of the be am. As discussed 
in §3.3 of Lopez-Canie go et al.l ()2007l ). this method for 
flux estimation through linear flltering is almost optimal 
in many circumstances. 

Note that the method adopted here to estimate the flux 
differs from the one used by the WMAP team. Assuming 
that they have followed similar procedures in intensity 
and polarization, they used a matched filter taking into 
account the non-Gaussian profile of the beam to detect 
point sources in the filtered full-sky maps, but their fluxes 
have been derived fltting the pixel intensities around the 
point source to a Gaussian profile plus a plane baseline 
(in the unfiltered image). 

We calculated the rms noise for each patch containing 
an input source. This value can be easily overestimated if 
border effects and strong fluctuations due to other point 
sources or small scale structure of the diffuse emissions 
in the patch are not removed or flltered out. In order 
to avoid this, first we find the 5% brightest pixels in the 
patch and flag them. Second, we flag pixels within a 
distance equal to 15 pixels from the border (see Fig. [3]). 
Finally, we select a shell around the source with inner 
radius equal to the FWHM of the beam (since polarized 
fluxes are never very high, this is enough for them to have 
decreased well below the noise level), and outer radius 
of 3 X FWHM. The rms noise, a, is then obtained as 
the square root of the variance of the pixels included in 
this shell, excluding flagged pixels (if any). From the 
distribution of the values of P in the unflagged pixels we 
calculate the levels exceeding those of 95.00%, 99.00%, 
99.90%, and 99.99% of pixels. These are the probabilities 
that signals at those levels are real rather than noise 
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Fig. 3. — These panels show an example of the flagging of the 
borders and the brightest 5% objects in a given patch to prevent 
contaminated pixels to be used in the calculation of the background 
distribution. The upper panel shows a filtered patch where the 
bright sources are identified. The lower panel shows, cross-hatched, 
the flagged area. 



TABLE 1 

Number of detections at different significance levels 



Significance 


23 GHz 


33 GHz 


41 GHz 


61 GHz 


95.00% 


138 


122 


93 


81 


99.00% 


53 


41 


30 


20 


99.90% 


20 


16 


11 


6 


99.99% 


18 


12 


9 


6 



fluctuations. An example is shown in Fig.O 
3. RESULTS 

In Table [1] we list the number of objects with P val- 
ues significant at more than 95.00%, 99.00%, 99.90%, 
and 99.99% confidence levels. There are 18, 12, 9, and 
6 detections with > 99.99% significance level at 23, 33, 
41, and 61 GHz, respectively. The 22 sources detected 
at such significance level at, at least, one frequency, are 
listed in Table [2] Several of them are well known bright 
extended objects including Cassiopeia A, Centaurus A, 
Cygnus A, Fornax A, Pictor A, Taurus A (Crab Nebula), 
and Virgo A (but, except for Cen A and For A, their 
sizes are smaller than the WMAP 23 GHz beam). Since 
our algorithm is not optimized for extended sources, the 



flux estimates for these objects are likely affected by sys- 
tematic errors much larger than the quoted statistical er- 
rors. In addition resolution effects can be seen in WMAP 
data, for example they likely account for the steep drop 
of the Crab polarized flux between 41 and 61 GHz. Such 
drop corresponds to the ratio of WMAP beam areas at 
these two frequencies. Five objects do not have plausible 
counterparts in low frequency radio surveys. They may 
therefore be exceptionally high peaks in the highly non- 
Gaussian fluctuation field, in temperature and polariza- 
tion, mostly due to Galactic synchrotron (with possible 
CMB and source confusion contributions). In fact, 3 of 
these 5 objects are at \b\ < 10° and a fourth one (no. 
346) is in the Ophiuchus Complex region. 

All the 5 obje cts detected in polarization by 
IWright et al.l (|2OO90 have P values significant at > 
99.99% levels at least at 1 frequency. As shown by Ta- 
ble [3] our flux estimates are in generally good agreement 
with those by IWright et al.l (2009) , in spite of the differ- 
ent techniques used. The main formal discrepancy con- 
cerns the 23 GHz polarized flux of 0322-3711 (Fornax 
A), a well known extended source for which we expect 
that photometric errors are mostly systematic, as noted 
above. 

To better assess the reliability of our estimates we have 
compared those sources detected at > 95.00% confidence 
levels with ground based measurements at nearby fre- 
quencies. Unfortunately there are only few samples that 
can be used for such a comparison. The AT20G Bright 
Source Sample (ESS; Massardi et al. 2008) covers the 
declination region S < —15° and is complete in total in- 
tensity down to S'20GHz — 0.5 Jy (except for Fornax A) 
and has simultaneous polarization measurements. Nine 
extended (i.e. with size larger than the 2.4 arcmin resolu- 
tion of ATCA measurements) objects have been foUowed- 
up in polarization, mosaicking a region large enough to 
evaluate reasonably well the integ rated polarized flux 
density (|Burke-Spolaor et al.l l2009f ) . Despite all efforts, 
for some very extended objects (like Cen A) the full ex- 
tent of the low-frequency radio structure could not be 
imaged and the total P could not be measured. Angu- 
lar resolution plays a key role in the comparison of the 
measurements of P obtained with different instrument or 
configurations for extended sources. The problems are, 
of course, amplified if observations at different frequen- 
cies are compared: the AT20G observations showed that 
the polarization degree may vary with frequency. The 
BSS comprises 218 detections of polarized flux density, 
of which 19 are above 100 mJy at 20 GHz. 73 objects 
of our 95.00% confidence level detections are in common 
with the BSS, but only 28 of them have a polarization 
detection in the AT20G BSS, with polarized flux density 
typically lower than those estimated with the techniques 
applied here, as detailed later. 

Ricci et al. (2004) have carried out 18 GHz po- 
larization ob s ervati ons of the Southern portion of the 
iKuehr et al] (|198lD sample, comprising sources with 
'S'sGHz = 1 Jy- Due to resolution effects, a source in the 
Centaurus A region appears with ffux density larger than 
1 Jy in our samp le but with less than 100 mJy in the 
iRicci et "all ()2004D sample. Variability may justify the 
disagreement for the blazar 3C 279 for which Ricci et al. 
measured Pisghz = 1-6 Jy, a factor of ~ 2.5 higher than 
estimated from the WMAP 23GHz map. 
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TABLE 2 

The polarized WMAP Point Sources (POWPS) sample. 



Object 


RA 


DEC 


GLON 


GLAT 


P (Jy) 


P (Jy) 


P (Jy) 


P (Jy) 


flags 




[h] 


[dcg] 


[dcg] 


[dcg] 


23 GHz 


33 GHz 


41 GHz 


61 GHz 





74 


(Fornax A) 


3, 


.372 


-37, 


.177 


240, 


.122 


-56, 


.766 


1.07 ± 0, 


,04 


0.87 


± 0, 


,06 


0.49 ± 0, 


,06 


0, 


.41 


± 


0. 


10 


1 


1 


1 


4 


126 


(Pictor A) 


5, 


,326 


-45, 


.743 


251, 


,548 


-34, 


,680 


0.38 ± 0, 


,05 


0.42 


± 0, 


,06 


0.41 ± 0, 


,07 






— 






1 


1 


3 





156 


(PKS0607-15) 


6, 


.166 


-15, 


.679 


222, 


.609 


-16, 


.109 


0.33 ± 0, 


,04 


0.36 


± 0, 


,07 
















1 


3 








200 


(PKS0829+04) 


8, 


.530 


4, 


.559 


220, 


.627 


24, 


,357 


— 






— 




0.56 it 0, 


,09 


1 


.10 


± 


0. 


16 








3 


1 


256 


(PKS1144-37) 


11, 


,787 


-38, 


.150 


289, 


,265 


23, 


,012 


— 




0.49 


± 0, 


,07 


— 








— 









1 








266 


(NC) 


12, 


.200 


-52, 


.630 


296, 


.882 


9, 


778 


0.73 ± 0, 


,04 


0.42 


± 0, 


,06 


0.44 ± 0, 


,07 












1 


1 


1 





272 


(3C273) 


iz. 


.4o0 


Z. 


,U44 


OQO 

zoy. 


.yo / 






l.U/ it U, 


■ UO 




it U, 


■ UD 


U. i z it U, 


HQ 

■ Uo 












1 


i 


1 


U 


273/274 (Virgo A) 


12, 


,509 


12, 


.350 


283, 


,597 


74, 


,433 


0.79 ± 0, 


,05 


0.71 


±0, 


,07 


0.62 ±0, 


,08 


0, 


,53 


± 


0. 


14 


1 


1 


1 


1 


280 


(3C279) 


12, 


.936 


-5, 


.762 


305, 


.107 


57, 


,090 


0.67 ±0, 


,05 


0.45 


±0, 


,08 


0.83 ±0, 


,09 


0, 


.66 


± 


0. 


17 


1 


3 


1 


4 


289 


(PKS1320-44) 


13, 


.381 


-44, 


.682 


308, 


.790 


17, 


,832 


1.67 ±0, 


,08 


1.13 


±0, 


12 


0.56 ±0, 


13 












1 


1 


3 





291 


(Centaurus A) 


13, 


,422 


-43, 


.025 


309, 


,483 


19, 


,416 


3.19 ± 0, 


,08 


2.30 


±0, 


11 


2.02 ±0, 


11 


1, 


,58 


± 


0. 


13 


1 


1 


1 


1 


295 


(3C286) 


13, 


.527 


30, 


.510 


56, 


.332 


80, 


,578 


0.33 ± 0, 


,04 


0.50 


±0, 


,06 






0, 


.63 


± 


0. 


11 


1 


1 





3 


337 


(PKS1546+02) 


15, 


.823 


2, 


.545 


10, 


.743 


40, 


,891 


0.33 ±0, 


,05 






















1 











346 


(NC) 


16, 


,345 


-25, 


.487 


351, 


,319 


17, 


150 


0.45 ±0, 


,06 






















1 











432 


(NC) 


20, 


.278 


45, 


.776 


82, 


.152 


5, 


,810 


1.16 ± 0, 


,07 


0.67 


±0, 


,07 


0.38 ±0, 


,07 












1 


1 


3 





437 


(GB6 J2038+51) 


20, 


,647 


51, 


.316 


88, 


.822 


6, 


,016 


0.62 ± 0, 


,07 


0.29 


±0, 


,07 


0.40 ± 0. 


,07 












1 


4 


3 





439 


(NC) 


20, 


.846 


29, 


.160 


72, 


.753 


-9, 


.460 


0.45 ±0, 


,05 


0.33 


±0, 


,06 


0.31 ±0, 


,07 


0, 


.48 


± 


0. 


11 


1 


3 


4 


4 


440 


(G93. 3+6.9) 


20, 


.871 


55, 


,404 


93, 


.307 


6, 


,956 


0.57 ±0, 


,06 


0.36 


±0, 


,07 


0.33 ±0, 


,06 












1 


4 


4 





473 


(NC) 


22, 


,323 


26, 


.438 


85, 


.400 


-25, 


,138 
















0, 


,95 


± 


0. 


13 











1 


Cyg 


nus A (3C405) 


19, 


,984 


40, 


.484 


75, 


.930 


5, 


,700 


0.49 ± 0, 


,07 


0.51 


±0, 


,07 


0.56 ±0, 


,06 


0, 


,48 


± 


0. 


12 


3 


2 


1 


4 


Taurus A (Crab) 


5, 


.583 


22, 


.369 


184, 


.310 


-5, 


.510 


24.7 ±0, 


,18 


20.2 


±0, 


,16 


16.0 ±0, 


,14 


6, 


.41 


± 


0. 


18 


1 


1 


1 


1 


Cas A (3C461) 


23, 


,408 


58, 


.835 


111, 


,872 


-2, 


,157 


0.91 ± 0, 


,05 


0.55 


±0, 


,06 


0.58 ±0, 


,06 


0, 


.22 


± 


0. 


11 


1 


1 


1 


1 



Note. — Sources detected at > 99.99% confidence level in at least one of the WMAP frequency channels. Column 1: sequential number 
in the NEWPS_5yr_5s catalogue and source name (NC means that the source has no plausible low radio frequency counterpart; see text); 
columns 2-5: equatorial (J2000) and Galactic coordinates of the source; columns 6-9: detected integrated polarized flux density and their 
errors at 23, 33, 41 and 61 GHz; columns 10-13: flags for significance (1 : > 99.99%, 2 : > 99.90% but < 99.99%, 3 : > 99.00% but 
< 99.90%, 4 : > 95.00% but < 99.00%, =< 95.00%), at 23, 33, 41 and 61 GHz, rospoctivoly. The NEWPS.5yr.5s catalogue lists 2 
sources (No. 273 and 274) close to the position of Virgo A. The present re-analysis has shown that they are actually the same source, 
coinciding with Virgo A with total fiux density of 18.4 ib 0.25 at 23 GHz. 



TABLE 3 

Comparison of polarized flux estimates in the present paper, with those of IWright et alI 

ipoogj) 



Object 23GHzw 23 GHzp 33GHzw 33GHzp ilGHz^ 41GHzp 

0322-3711 (For A) 1.57 ±0.07 1.07 ±0.04 1.17 ±0.15 0.87 ± 0.06 0.85 ± 0.23 0.49 ± 0.06 

0519-4546 (Pic A) 0.39 ± 0.06 0.38 ± 0.05 0.45 ±0.10 0.42 ± 0.06 - 0.41 ± 0.06 

1229+0203 (3C273) 0.98 ± 0.06 1.07 ±0.05 0.81 ±0.11 0.71 ±0.06 0.80 ±0.13 0.72 ± 0.08 

1230+1223 (Vir A) 0.75 + 0.08 0.79 ±0.05 0.74 ±0.11 0.71 ±0.07 0.50 ±0.11 0.62 ± 0.08 

1256-0547 (3C279) 0.62 ± 0.07 0.67 ±0.05 0.55 ±0.12 0.45 ± 0.08 0.69 ±0.15 0.83 ± 0.09 



Note. — The subscripts 'W and 'p' designate the columns containing the results bv lWright et all (120091) and 
of the present paper, respectively. 



Both the BSS and the Ricci et al. samples cover the 
Southern hemisphere. Among the Northern 78 VLA po- 
larization calibrators listed in the compilation updated 
by S. Meyer'^ 31 are among our 95.00% significance level 
detections. 13 of them have polarized flux density (aver- 
aged over the epoch 2002-2006) above 100 mJy in the K 
band. 

As mentioned above, the comparison with ground 
based measurements is complicated by resolution effects 
for the extended (generally steep-spectrum) sources and 
by the strong variability of the flat-spectrum ones. Nev- 
ertheless, Fig. m shows a reasonably good consistency for 
P > 400 mJy at 23 GHz. Below 400 mJy, the values of P 
are clearly dominated by the contribution of positive po- 
larization fluctuations at the source positions (Eddington 
bias). 

3 www.aoc.nrao.edu/ smyers/calibration/master.shtml 



3.1. Additional tests and simulations 

Unfortunately the number of objects with polarized 
flux above 400 mJy and ground based polarization mea- 
surement at frequencies close to WIVIAP ones is small, 
so that our comparison has a poor statistics. To better 
assess the reliability of our flux estimates it is thus neces- 
sary to resort to simulations. We have selected a sample 
of 738 positions with |6| > 5° and far from each other 
more than ARg (see the first paragraph of ii l2.3p and more 
than 2Rs away from each object in the NEWPS_5yr_3s 
catalogue. These 'blank' positions (free of sources 
brighter than 3(t in total intensity) constitute our con- 
trol fields. We have then chosen 10 values of P, ranging 
from 0.2 to 2 Jy, with a constant step in logP, i.e. P = 
[0.2, 0.26, 0.33, 0.43, 0.56, 0.719, 0.928, 1.20, 1.55, 2.0] . For 
each value of P we have injected a source with that po- 
larized flux density in the projected Q and U patches 
centered on 100 randomly chosen control field positions. 
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Fig. 4. — Comparison of polarized flux density estimates from the WMAP 23 GHz map (present paper) and ground-based observations: 
AT20G BSS (Massardi et al. 2008; squares); Ricci et al. (2004) (diamonds); VLA calibrators (triangles). Filled symbols correspond to 
sources detected on the WMAP map with more than 99.99% significance level. In case of sources with multiple ground-based observations, 
we have chosen those with the resolution closest to WMAP's. The solid line corresponds to equal values on the two axes. The two highest 
filled symbols on the left of the solid line are examples of the effect of the much higher resolution of ground based measurements, compared 
to WMAP's; the filled diamond is Cen A (resolved also by WMAP), and the filled triangle is Virgo A (unresolved by WMAP). The highest 
open triangle corresponds to a source close to the Galactic plane, whose polarized flux estimated from the WMAP maps includes a dominant 
Galactic contribution. The outlier on t he right of the solid line is the highly variable blazar 3C279. The asterisks with error bars show the 
results of simulations described inS IS.ll 



This was done randomly selecting Q between —P and 

P and setting U = ±^/P^ - Q'^, the sign of U being 
again chosen at random. The inserted source was con- 
volved with the WMAP symmetrized beam at 23 GHz. 
We avoided using the same patch more than once, except 
when all the control fields were already used. In any case 
the same patch was never used twice for the same value 
of P. 

Figure [5] shows the percentage error in the values of 
P of the simulated sources recovered with our filtering 
and flux estimation process with respect to the input val- 
ues. This comparison confirms the reliability of recovered 
fluxes for Pinput ^ 400 mJy, while the Eddington bias 

becomes increasingly important below this value. For 
^input = 400 mJy the recovered fluxes are, on average, 
overestimated by less than 10%. 

3.2. Median Polarization Degree 

The median polarization degree at 23 GHz of the 11 
sources detected with a confldence level > 99.99%, with 
low radio frequency counterparts, and with P > 400 mJy 
is ~ 7.5%. This value, however, is not representative of 
the mean polarization level of sources in a complete sam- 
ple for two reasons. First, the sample is obviously biased 
towards sources with the highest polarization degrees. 



10 



5 - 




-5- 



0.5 1 1.5 2 2.5 

P 

Fig. 5. — Results of the simulations for logarithmically spaced 
values of P in the range 0.2-2 Jy. In this panel we show the percent- 
age error (100( 

-^input -^recovered )/^input) in the recovered value 

of P. 



Second, as mentioned above, our photometry is not op- 
timal for extended sources which make up a substantial 
fraction of the sample, so that their estimated polariza- 
tion degrees are highly uncertain. 
An unbiased estimate of the mean polarization degree 
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Fig. 6. — Distributions, normalized to unity, of P signals in the di- 
rection of the 23 GHz sources brighter than 5 Jy (13 objects; dashed 
histogram) compared with those for the control fields (solid). The 
median P values are of 0.27 Jy for sources and of 0.16 Jy for control 
fields. According to the Kolmogorov-Smirnov test, the probability 
that the two samples are drawn from the same parent distribution 
is of only 7.3 • 10"^ 



of sources at WMAP frequencies can be obtained from a 
comparison of the distribution of P values for a suitably 
chosen complete subsample with that of control fields. 
The completeness limit, in total flux, cannot be too faint, 
otherwise the mean polarization level is too low to be 
detectable. On the other hand, if the flux limit is too 
high, the number of sources is too small for a meaningful 
statistical inference. The optimal flux limits turn out 
to be of 5, 4, 3, and 4 Jy at 23, 33, 41, and 61 GHz, 
respectively. We find a highly significant detection of 
polarized flux density only at 23 GHz (see Fig. [6]): the 
probability that the distribution of P values of sources is 
drawn from the same parent distribution as control flelds 
is 7.3 • 10~^. The median polarization degree, H, can be 
estimated as: 



H 



P 



mod, sources 



P 



incd.coiitrolficlds 



5med 



(3) 



sources 



where ^med, sources is the median flux of sources and we 
have neglected the median flux of control flelds, which 
is close to zero. The result at 23 GHz is H23GHZ = 
1.7 ± 1.1%, consistent with the median polarizat ion de- 
gree f or the AT20G BSS, Hsoghz - 2.5% (MassardieLalJ 
12008( 1 ■ Note that the lower resolution of WMAP, com- 
pared to AT20G, observations make them more liable to 
beam depolarization (due to chaotic components of the 
magnetic fleld within the unresolved region) in the case 
of extended sources. 

The probabilities that the distributions of P values of 
sources are drawn from the same parent distribution as 
control flelds are of 0.02, 0.012, and 0.065 at 33, 41, 
and 61 GHz, respectively. The derived values of H are 
H33GHZ = 0.91 ± 0.83%, H41GHZ = 0.68 ± 1.0%, and 
HeiGHz - 1.3 ±1.8%. 

3.3. Predictions for Planck 

The Planck satellite is expected to measure the po- 
larization of the sources with greater sensitivity than 



WMAP, and therefore to detect more sources, down to 
fainter polarized flux limits. Building on the results of 
our analysis of WMAP maps, we have carried out simula- 
tions to estimate the minimum polarized flux detectable 
and reliably measurable by Planck. We have adopted 
the nominal mean instrumental noise levels expected af- 
ter two complete sky surveys (1 yr)'*. 

Assuming that the Q and U images are dominated by 
instrumental white noise, and adopting a pixel size of 
6.87 arcmin and idealized matched fllters for the nomi- 
nal beam sizes, we have computed the ratios between the 
(7p levels for the 30-100 GHz Planck frequency chan- 
nels and those for the closest WMAP frequency chan- 
nels. We flnd Crp,W?MAP5yr/crp,Plancklyr = 2.2, 1.6, 2.0, 

and 6.8 at about 30, 44, 70, and 100 GHz, respectively. 
Note that the higher Planck sensitivity is partly compen- 
sated by the longer WMAP exposure time. An extension 
of the Planck mission for one more year would decrease 
o-ppianck by a factor %/2. 

The above calculations take into account only instru- 
mental noise. To investigate the effect also of polariza- 
tion fluctuations due to diffuse foregrounds and to the 
CMB we have performed simulations analogous to those 
we did for WMAP at 23 GHz for the Planck 30 GHz 
c hannel, using the Planck Sky Model simulation software 
(Del abrouille et al.|[2009l ). 

In simulated Q and U maps containing polarized dif- 
fuse foregrounds and the CMB we have injected sources 
with P in the range 0.05-0.6 Jy. Our algorithm recov- 
ered sources with P down to 300 mJy with systematic 
offsets due to the Eddington bias of few percent or less. 

This result already highlights the difflculty of flnding 
suitable polarization calibrators for Planck. Not many 
point-like sources have P substantially larger than 300 
mJy at 30 GHz. Combining the results by Ricci et al. 
(2004), of the AT20G BSS (Massardi et al. 2008) and 
the VLA calibrator observations, we have found 11 ex- 
tragalactic sources above this limit. The count may be 
incomplete, particularly in the Northern hemisphere, but 
it is unlikely that many bright sources have been missed. 
Our analysis of 33 GHz WMAP maps, including bright 
Galactic sources, has detected 12 sources at a confldence 
level > 99.99%; polarized flux densities for most of them, 
however, are overestimated because of the Eddington 
bias. Also, some of the brightest sources are extended 
(compared to the Planck beam) and others (3C273 and 
3C279) are highly variable. 

4. CONCLUSIONS 

We have applied to the WMAP 5-yr polarization maps 
a new source detection technique, called "flltcrcd fusion" 
(Argiieso et al. 2009), taking into account the real sym- 
metrized beam profiles. The technique was applied at 
the positions of WMAP sources detected at > 5a by 
Massardi et al. (2009), plus 3 extended sources (Cygnus 
A, Taurus A, Gas A) known to be bright and polarized. 
22 sources were detected at a confidence level > 99.99% 
(that would correspond to ~ 3.72cr for a one-tailed Gaus- 
sian statistics) in at least one WMAP channel. Five of 
them, however, have no plausible counterparts in low ra- 
dio frequency catalogs and may therefore be just high 
intensity peaks of the fluctuation fleld. Nevertheless, 

4 www.rssd.esa.int/index.php?project=Planck 
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this is a substantial improvement compa red to the 5 
source polarization measurements listed by I Wright et alj 
((2001 and our results for these 5 sources are in gener- 
ally good agreement with theirs. There are several reasons 
that could explain such an improvement. First, WMAP 
used a mask to produce the point source catalogs based 
on the KQ75 mask plus the Magellanic cloud regions that 
exclude 11 objects that we have studied. Second, WMAP 
did not consider objects that are not associated to known 
sources at low frequencies. In our case, we do a non-blind 
search at the positions of the NEWPS-5yr catalog that 
includes a handful of such objects, some of which show 
a significant polarized flux (marked as NC in Table [5]) , 
maybe due to Galactic emission in polarization. Taking 
all these into account, and considering the efficiency of 
our method as compared to WMAP, one could conclude 
that we are detecting between 30 — 40% more sources, 
which is a similar improveme nt to the one obtaine d in 
our analysis in total intensity (jMassardi et al.|[2008[ ). 

A comparison of our polarized flux, P, estimates at 23 
GHz (Fig. m where sources detected at > 99.99% confi- 
dence level are represented by filled symbols) with high 
signal-to-noise ground based measurements at nearby 
frequencies, highlights the complications due to differ- 
ent angular resolutions in the case of extended sources 
(Gen A and Vir A) and of strong variability (3G279). 
For sources not affected by these problems, the agree- 
ment is quite good. The latter sources, however, are too 
few to allow a firm conclusion on the reliability of our flux 
estimate. We have therefore resorted to simulations, in- 
jecting fake sources of known polarized flux density in the 
WMAP 23 GHz Q and U maps. The simulations showed 
that our approach is reliable for P23GHZ > 400 mJy, while 
the Eddington bias becomes increasingly large at fainter 
fluxes (it is less than 10% at 400 mJy). Figure [4] also 
shows that estimates of P for sources detected at a con- 
fidence level < 99.99% are badly affected by the Edding- 
ton bias, and therefore unreliable. This was expected 
since, for example, a 99.90% confidence level corresponds 
to 3.1f7 in the case of a Gaussian distribution of fiuctu- 
ations. And, as shown by Hogg & Turner (1998), fiux 
estimates at this confidence level are practically useless. 

Analogous simulations using the mean noise lev- 
els expected after 1 year of Planck observations have 
shown that a few percent Eddington bias is reached for 
P30GHZ — 300 mJy. The modest difference compared to 
WMAP 23 GHz is due to the fact that the longer WMAP 
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exposures (5 yrs vs. 1 yr) compensate for the lower sen- 
sitivity. In the case of an extension of the Planck mission 
for a second year, the gain over WMAP will increase by 
a factor of Based on the predictions by Taylor et 
al. (2007; their Fig. 7), this would imply an increase of 
detections in polarization by a factor ^ 3 compared to 

WMAP. 

Sources detected on WMAP polarization maps have, 
on average, exceptionally high polarization degrees be- 
cause only a very small fraction of sources have been 
detected and they likely populate the tail of the distri- 
bution of P values. Estimates of, or upper limits to, 
the mean polarization degrees, H, of bright sources at 
23, 33, 41, and 61 GHz have been obtained comparing 
the distributions of polarized flux densities in the direc- 
tions of complete source samples, limited in total flux, 
with those in the directions of control fields, devoid of 
bright sources. The derived values, H23GHZ = 1.7±1.1%, 
H33GHZ = 0.91 ± 0.83%, H41GHZ = 0.68 ± 1.0% and 
HgiGHz = 1.3 ± 1.8%, are consistent with, but somewhat 
on the low-side of, the median polarization deg ree for 
the AT20G BSS, H20GHZ ^ 2.5% (IMassardi et al.i r20081. 
Somewhat lower values of H are expected at the WMAP 
resolution for extended sources, due to beam depolariza- 
tion. 

Finally, the detected sources may be candidate calibra- 
tors for other high sensitivity GMB experiments, like the 
ground-based QUIJOTE experiment (Rubiiio-Martin et 
al. 2008) or the Planck satellite mission. 
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